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March 1, 2011 

Abstract Two aspects of bright matter-wave solitons 
in weak external potentials are discussed. First, we briefly 
review recent results on the Anderson localization of an 
entire soliton in disordered potentials [I] , as a paradig- 
matic showcase of genuine quantum dynamics beyond 
simple perturbation theory. Second, we calculate the 
linear response of the mean-field soliton shape to a 
weak, but otherwise arbitrary external potential, with 
a detailed application to lattice potentials. 

PACS 03.75.Lm Tunneling, Josephson effect, Bose- 
Einstein condensates in periodic potentials, solitons, 
vortices, and topological excitations • 05.60.Gg 
Quantum transport 



1 Introduction 

Recently, the cold-atom community has shown renewed 
interest in soliton dynamics, sparked by the experimen- 
tal observation of cold-atom solitons in quasi-oncdimcn- 
sional Bose-Einstein condensates with attractive con- 
tact interaction [HE]. Notably, it has been emphasized 
that the soliton's center of mass is a collective degree of 
freedom whose dynamics can show genuine quantum ef- 
fects. In this vein, Weiss and Castin [3] have calculated 
the scattering amplitude of a soliton by a potential bar- 
rier, which results in a superposition of classically dis- 
tinct quantum states, namely the soliton being either 
transmitted or reflected. Similarly, Lewenstein and Mal- 
omed 5i have proposed to generate entanglement by 
the controlled collision of quantum solitons. 

C. A. Miiller, Centre for Quantum Technologies, National Uni- 
versity of Singapore, Singapore 117543, Singapore 
E-mail: cord.mucllcr@nus.cdu.sg 



A well-known paradigm of genuine quantum dy- 
namics is disorder-induced Anderson localization [51[7], 
which has been studied for solitons in different settings 
some time ago [5], and also been observed rather re- 
cently with ultracold, noninteracting matter waves 
fTUlfTTUTS] (see also [H]). Motivated by these experimen- 
tal advances, we have investigated the quantum dynam- 
ics of matter-wave solitons in spatially correlated dis- 
order potentials [1] . The first part of the present paper 
reviews briefly the derivation of an effective Hamilto- 
nian for the center of mass and the resulting localization 
exponent in an optical speckle potential. 

Naturally, there is more to solitons than just their 
center-of-mass dynamics. Whenever inhomogeneous force 
fields act on a compound object, the latter responds 
by adapting its internal configuration as well. For weak 
forces, the response will be linear and thus described by 
a (linear) susceptibility. The second part of this paper 
investigates in more detail how the soliton's ground- 
state shape changes under the influence of a weak ex- 
ternal potential. We calculate the linear compressibility 
in general, and then focus on the simple, yet interesting 
case of a lattice potential, including a detailed compar- 
ison of analytical results to numerical data. This allows 
us finally to derive quantitative criteria for the external 
perturbation of the shape to be weak. 



2 Setting the stage 

We describe a weakly interacting Bose-Einstein conden- 
sate (BEC) in a quasi-onedimensional wave guide by its 
mean-field amplitude (f>(z). The Gross-Pit aevskii (GP) 
free energy functional, with a given chemical potential 



fi and in a homogeneous wave guide, reads 

Eofaf] - |dz |i|9 2 0| 2 + ||^| 4 - mI^I 2 } ■ (1) 

We will use units such that K = m = 1 throughout 
the following, g = 2uj_a is the effective interaction 
constant for a quasi-onedimensional condensate with s- 
wave scattering length a and transverse harmonic trap- 
ping frequency uj±. The cases of repulsive (a > 0) and 
attractive (a < 0) interaction correspond to g > and 
g < 0, respectively. 

Minimizing the free energy Eq yields the ground 
state 4>o(z). The chemical potential [i thus determines 
the total number of particles A*o = J dzno(z), with 
n (z) = \4> (z)\ 2 being the condensate density. Let us 
choose periodic boundary conditions, as for a toroidal 
wave guide with circumference L. In the case of repul- 
sive interaction g > 0, both kinetic energy and interac- 
tion energy are mimimized by spreading the density ho- 
mogeneously over the entire available length: no = fi/g 
and No = L[ifg with /j, > 0. The corresponding ground- 
state wave function O = e j6, ° \fjjjg plays the role of 
a BEC order parameter, featuring an arbitrary global 
phase 8 that spontaneously breaks the U(l) gauge in- 
variance of (fil. For the ground state of a single Bose- 
Einstein condensate, we can set 8q = 0. 

An attractive interaction g < rather favors a state 
where atoms are clustered together. As shown by Kana- 
moto et al. [14) . for large enough system size or chemical 
potential (i.e. number of particles), Ll/i] 1 / 2 ^> 1, the 
energetically preferred state is a soliton, 
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(2) 



with a hyperbolic secant envelope, decaying over the 
length scale £ = |2/x| -1 / 2 known as the condensate 
healing length. The number of particles is now No — 
2|2/i| 1 / 2 /|g|. Conversely, in a canonical setting with fixed 
number of particles N, the chemical potential settles to 
Mo = -V(2£ 2 ) = -g 2 N 2 /8 < 0. 

Note that the solution |2]) spontaneously breaks not 
only the gauge invariance with a phase 9q (that we 
take to be zero), but also the translational invariance 
of the energy functional ([I]). Therefore, the center of 
mass zq ='■ q emerges as a dynamical degree of freedom 
on its own. In the homogeneous situation described by 
(II]), all classical solutions q(i) = qo + vt with constant 
velocity v are admissible by Galilean invariance. In a 
given frame of reference, the classical configuration with 
minimum kinetic energy is a soliton resting at q = qo- 
Quantum mechanically, however, q is rather distributed 
with its ground-state wave function &o(q) that obeys 



the free Schrodinger equation under the given bound- 
ary conditions. For periodic boundary conditions, this 
center-of-mass ground state is the completely dclocal- 
ized plane wave with momentum p q — 0, i.e. the con- 
stant &o{q) = L~ x / 2 . Contrary to the point of view that 
the "wave properties of solitons manifest themselves in 
radiation emitted due to scattering by impurities" [51 
p. 53], we take a more general approach along the funda- 
mental principles of quantum mechanics, namely that 
wave-particle duality is a general feature of all degrees 
of freedom, to be revealed under appropriate experi- 
mental circumstances. 

It may seem peculiar to speak of quantum dynamics 
for the center of mass of a whole collection of atoms, 
after starting out from a mean-field description for the 
entire condensate. But it is indeed very common to 
separate the center of mass from internal variables in 
interacting systems and moreover exact for two-body 
forces that only depend on the relative distance between 
the microscopic constituents. The situation of a soliton, 
composed of individual atoms held together by attrac- 
tive contact interactions, is therefore analogous to the 
situation of, say, an alkali atom composed of a nucleus 
and electrons held together by the Coulomb force. And 
just as one may study the quantum dynamics of en- 
tire atoms — not to mention fullerenes or biomolecules 
|15| — one may indeed equally well study the quantum 
dynamics of entire solitons as far as their center of mass 
is concerned. 

Consider now the mean-field energy functional in 
presence of an external potential V(z): 



E = / dz 



,\d z ^ 



A^ + [V{z)-nM 



(3) 



We assume that V(z) is a small perturbation on the en- 
ergy scale set by /i. To lowest order in V/|/x|, therefore, 
the shape of the original soliton will remain unchanged. 
In contrast, the external potential can never be consid- 
ered a small perturbation for the center-of-mass dynam- 
ics because a finite V is never small compared to the 
homogeneous case V = 0. 

At this point, the paper bifurcates. In the following 
section, we focus on the center of mass. We highlight the 
quantumness of its dynamics in a disordered potential 
by discussing the Anderson localization length, as first 
derived in PQ. As a complement to this work, we will 
study in sections [4] and [5] below how the soliton's shape 
is modified by the presence of a weak external potential 
in the mean-field ground state. Section [6] concludes. 



3 Anderson localization of a soliton 

3.1 Effective Hamiltonian 

Assuming a fixed soliton shape, the center of mass q 
of TV particles can be described as a collective variable 
with the ansatz cj)(z; q,p q ) = e lPqZ / N (J)q(z — q). Here, the 
conjugate momentum p q appears in the phase together 
with a factor N~ x because 4>o, eq. (pL is normalized 
to N = 2|2/i| 1 / 2 /|<?|. Inserting this collective-variable 
ansatz in the energy functional © and integrating over 
z yields the effective Hamiltonian 



L-ifdz 



-ikz 



V(z). Statistical homogeneity then trans- 



q ~ 2N 



&z\Mz-q)\ 2 V(z). 



(4) 



This Hamiltonian describes a particle with mass N evolv- 
ing in a potential V(q) = V * no{q) that is the convolu- 
tion of the bare potential with the soliton envelope. 

If the bare potential varies only very slowly over 
one healing length £ = 2/(N\g\), then the soliton feels 
the sum of forces on its constituents, V(q) — NV(q). If 
on the other hand the potential varies rather rapidly, 
then the convolution by the soliton density washes out 
all details on scales smaller than £, and the effective 
potential is strongly reduced. This is easily illustrated 
with a lattice potential V(z) — Vo cos(kz), for which the 
effective potential is essentially the Fourier transform of 
the soliton density \(f>o( z — z a)\ 2 of (pi): 



V(q) = N 



?rfc£/2 



sinh(7rfc£/2) 



Vocos(kq). 



(5) 



Indeed, V{q) -> NV(q) as fc£ -> 0, and for fc£ > 1, the 
amplitude Vq ~ ke^^^l 2 becomes exponentially small. 



3.2 Correlated disorder 

A disorder potential is a random process V(z) char- 
acterized by its moments V(z), V(z\)V{zz), etc. Sta- 
tistically homogeneous disorder is translation invariant 
after averaging, with V(z) — const., V(z + z )V(z ) — 
V(z)V(0), etc. Without loss of generality, one can al- 
ways set V(z) — by redefining the zero of energy or 
shifting n i->- \x + V(z) in pj. Thus, the most basic in- 
formation about the disorder potential is its pair corre- 
lator V(z)V(0) = V$C(z/a) where V£ := V{z) 2 is the 
variance characterizing the overall strength of disorder. 
The spatial correlation function C(z/a) decreases from 
C(0) = 1 to zero over a characteristic length scale a, 
the correlation length of the disorder. 

Equivalently, a disordered potential can be seen as a 
random superposition of plane-wave components Vu — 



lates into conservation of total momentum under aver- 
aging: VfeVfe' = L~ 1 5fe j _fc'y o 2 P(/c), where the so-called 
power spectrum P(k) is the Fourier transform of the 
real-space correlator C{zja). 

For the present case of cold-atom dynamics, we con- 
sider in detail optical speckle potentials [16] for which 
the laws of optics result in a remarkably simple corre- 
lation: C{z/(j) = [sinc(z/er)] 2 or 



P{k) = ira{l - ||H)e(l - h ka \)- 



(6) 



The correlation length a is determined by the wave 
length of the laser light and the geometric aperture of 
the imaging system and can be as short as a — 0.26 /iin 
[9]. The Heaviside distribution (9(.) in (pi) excludes all 
wave vectors with modulus larger than 2/ct, as required 
by the limit of optical resolution. 

The effective potential felt by the soliton's center 
of mass is the convolution of the bare potential by the 
soliton density. By virtue of the convolution theorem, 
the Fourier components of the effective potential are 
therefore the product of the bare components times the 
Fourier components of the density, which already ap- 
peared in ([5]): 



V k = N 



rfc£/2 



sinh(7rfc£/2) 



V k . 



(7) 



The statistical properties of the potential affecting 
the soliton's center of mass are therefore completely de- 
termined and readily expressed in Fourier components. 
For example, the effective power spectrum reads 



P(k) = JV S 



(^£/2) 2 
sinh(7r/s£/2) : 



; P(k). 



(8) 



If the potential is very smooth on the soliton scale £, 
i.e. has a correlation length a ^$> £, the power spec- 
trum P(k) goes to zero faster than the soliton Fourier 
envelope, and P(k) ~ N 2 P(k). Conversely, if the bare 
potential is varying very rapidly, a <C £, the bare spec- 
trum can be approximated in the range ka ^C 1 by its 
delta-correlation limit P(0). Then, the healing length £ 
takes over as the new correlation length with an expo- 
nential decay of potential fluctuations as 

P(jfc) w {NTTk£) 2 P{0)e-^ (9) 

for fc£ > 1. 



3.3 Anderson localization exponent 

The Hamiltonian of the free soliton, p 2 /2N, has the 
plane- wave eigenfunctions &k(<l) ^ e lkq . When a small 




Fig. 1 Log-linear plot of the soliton localization exponent 7(fc) 
versus its center-of-mass wave vector k, in units of the soliton 
width £. Circles are numerical results obtained by diagonalizing 
the Hamiltonian Bl and solid lines are the result of a transfer- 
matrix calculation [T1|20| . for a red- and blue-detuned speckle 
potential with Vo = ±8 • 10 — 5 |/i|, respectively, and a correlation 
length a = 0.28£ ~ 0.26 fj,ta. The Born approximation, Eq. l |ll| , 
is shown as a dashed line. The overall exponential decrease for 
kt; ^> 1 and ka -C 1 is clearly visible. 



disorder potential is switched on, these extended plane- 
wave states become exponentially localized. Mathemat- 
ically rigorous theorems assure that in ID, the loga- 
rithmically averaged eigenfunctions decay for large dis- 
tances q from their origin like 



1 



lim log|<f(g)| = -- 7 (%. 

q—too Z 



(10) 



The corresponding localization exponent, or inverse lo- 
calization length, j(k) can be calculated perturbatively 
(for generic values of k, i.e., away from the band cen- 
ter k = or other singular points) in powers of the 
strength of the disorder potential. To second order in 
Vo, in the so-called Born approximation, the localiza- 
tion exponent reads j(k) = (/c 2 V 2 /4£|)P(2fc), with 
E k = h 2 k 2 /2m the free kinetic energy [T71I18U19J . The 
potential correlator is evaluated at momentum 2k since 
it is the elementary backscattering process k — > — k 
that is eventually responsible for Anderson localization 
in ID. 

For the soliton, we have Ek = k 2 /(2N) and scatter- 
ing by the effective potential Vfe, such that the localiza- 
tion exponent is predicted to be 



7(A) 



N 2 V 2 
k 2 



P(2k) 



N A V 2 (irkQ 2 
k 2 



; P(2k). 



(11) 



sinh(7rfc£) 2 

Figure [T] shows this prediction, together with numeri- 
cal data, obtained both by exact diagonalization of the 
Hamiltonian (El) and a transfer matrix approach, re- 
spectively [TII2T)]. We have chosen realistic experimen- 
tal parameters: N — 100 7 Li atoms with scattering 



length a = — 3nm in a transverse trap with uj± = 
2irx5 kHz form a soliton of size £ ~ 100 /im/TV w 1 /im. 
We consider an optical speckle potential with ampli- 
tude V = ±8 • 10~ 5 |^|, i.e. both the red-detuned case 
with Vo < an d the blue-detuned case with V > 0. 
Since the speckle potential has a non-Gaussian, skewed 
distribution, the full localization exponent depends on 
the absolute sign of Vo, an effect that the lowest-order 
Born approximation 0(Vq 2 ) cannot capture. However, 
the overall exponential decrease for fc£ 3> 1 is correctly 
predicted. 

These results show that the localization length, 7 -1 , 
of not-too-fast solitons with fc£ w 1 is in the sub-mm 
range that is measurable in current experiments on lo- 
calization of non-interacting matter waves [5] . Further 
results can be found in jT] , such as the expected density 
distribution of the final soliton position, depending on 
initial trapping conditions. 

But what about the emission of radiation due to 
scattering by impurities, as discussed at length by Gre- 
deskul and Kivshar [Sj and others? It is important to 
realize that we specifically focus on smooth, spatially 
correlated potentials that are expected to provoke con- 
siderably less excitation than isolated, (5-likc impurities. 
Nonetheless, it is clearly important to know the precise 
effects of the external disorder on the soliton shape. 
Some light on this issue will be shed in the following 
sections. 



4 Linear response for the mean-field ground 
state 

The previous results on the disorder-induced localiza- 
tion of a soliton have been derived with the effective 
Hamiltonian Q as the only ingredient. There, the soli- 
ton shape was assumed to be completely unaffected by 
the external potential. Although \V\ <C |/n| is probably 
a sufficient condition for this to be valid [20], as such 
it cannot be a satisfying, namely necessary criterion. 
Indeed, a potential with a large amplitude V » \fi\ but 
fluctuating only on very long length scales hardly de- 
formes the soliton if its variation over the soliton size, 
£dV/dz, is very small. Conversely, the amplitude of an 
interacting condensate cannot follow very rapid poten- 
tial oscillations with a wavelength much smaller than £ 



since this would cost too much kinetic energy in (12). 



In this case, only strongly smoothed fluctuations are 
expected to appear, just as in the repulsive case [2T] , 
In order to sharpen the picture, we therefore set out to 
compute the mean-field soliton deformation as function 
of the potential's amplitude and wave vector. 

A small deformation of the ground state can be de- 
scribed by an expansion over a complete set of elemen- 



tary excitations, such as described by Bogoliubov the- 
ory. In general, the Bogoliubov theory of condensate 
excitations in inhomogeneous potentials can be formu- 
lated as a saddle-point expansion of the mean-field en- 
ergy functional [22,23] and comprises two important 
steps: In a first step, the deformed ground state is de- 
termined as a functional of the external potential. In a 
second step, the quadratic excitations around this de- 
formed ground state are determined. In order to arrive 
at meaningful results with all effects to a given order 
of V taken properly into account, the first step is a 
vital prerequisite for carrying out the second. In the 
present contribution, we shall content ourselves with 
taking step 1, the consistent calculation of the mean- 
field ground state deformation caused by a weak, but 
otherwise arbitrary, external potential V(z). 

For the repulsive case with its homogeneous ground- 
state density no = fJ-/g, such a linear response to an ex- 
ternal inhomogeneous potential can be computed rather 
easily, taking advantage of Fourier decomposition; for 
applications to disordered potentials see, e.g., the work 
of Giorgini et al. [H] and Sanchez- Palencia [3T]. For a 
soliton with its inhomogeneous ground state density, a 
very similar linear response scheme, technically slightly 
more demanding but perhaps also more interesting, is 
presented in the following. 



4.1 Density response 

Taking our cue from Giorgini et al., we start with the 
GP functional (IT]) in density-phase representation <j){z) — 



y/n{z)e 



E 



i6(z). 



Az 



(a* 



n(d z ef 



9 2 

-n -[in 



(12) 



Since the perturbation J dV (z)n(z) couples to the den- 
sity, the deformed ground-state density reads n(z) = 
no(z) + 5n(z), where the shift within linear response is 
given by 



5n(z) 



dz'x{z,z')V(z')+0(V 2 ). 



(13) 



The density-density susceptibility (essentially the com- 
pressibility) is defined via its functional inverse, 



X^M 



S 2 E a 



Sn(z)Sn(z') 



2n (z) 



d z d z > 



n (z) 



3<? 



(14) 

S(z-z'). 
(15) 



Unlike in the repulsive case with no constant where 
the inverse is easily computed in Fourier modes [24] . 
here we have an inhomogeneous density Uq(z) and thus 
X -1 difficult to invert. We bypass this difficulty with 
a second approach in the upcoming section and come 
back to the compressibility in section [4~3| below. 



4.2 Wave-function response 

The mean- field ground-state wave function (f>(z), which 
can be taken real, satisfies 5E/S<fi* = 0, known as the 
stationary GP equation 



1 



d 2 + g(f>(z) 2 -n+ V(z) 



(z)=0. 



(16) 



Following Sanchez-Palencia |21J . we develop 4>(z) — 4>q{z- 
zq) + S(j)(z — Zq) around the ground-state solution (pi) 
by treating 5(f> as a small quantity of order V/|/x|. Lin- 
earizing gives 



1 



d 2 + 3gn (z - zq) - \x 



(z-z ) = -V(z)4>o(z-Zq). 
(17) 



Measuring distances in units of £ around z , i.e. writing 
z — zq = £ x and dividing by 2|u| = £~ 2 , one finds the 
linear equation 



— d 2 T -3scc1i(.t) 2 + - 
2 2 



dtp(x) = — tp (x) 



2|/x| 



(18) 



for the dimensionless shift Sip(x) = \/^S(j)(^x). This 
equation is of the form 



Ha 



1 



5ip(x) = W(x) 



(19) 



with Hq = p 2 /2 — 3sech(x) 2 the Hamiltonian of the 
sech -potential well, the well-known Poschl- Teller po- 
tential [25] . On the right-hand side, one has an external 
source term 



W{x) = -tp {x) 



V(zq + &) 
2|Ml 



(20) 



Computing 8(p then only requires to invert Hq. For- 
tunately for us, this is very simple to do because all 
potentials of the form —~i/(y + l)sech(x) 2 with inte- 
ger v are supersymmetric partners of the free-particle 
case v = 0, and their eigenstates and eigenenergies are 
perfectly known; for a brief and pedagogical introduc- 
tion to these issues, see 1:26 . The case v = 2 of interest 



here is treated in detail by Lekner [27]. Hq admits two 
bound states, 



ipi(x) 



sech(x) 2 , 



sech x tanh x 



(21) 
(22) 



with eigenenergies Eq = — 2 and E\ = — ^ (in units of 
2|/i|), respectively, as well as scattering eigenstates 



4>k{x) 



„ikx 



[27T]V2 



k 2 - 2 + 3 sech(x) 2 + iik tanh x 

[(l + fc 2 )(4 + /c 2 )] 1 /2 



(23) 



with free kinetic eigenenergy E k = fc 2 /2, k £ R. These 
eigenmodes have appeared repeatedly in one form or 
another in the soliton literature (cf. a recent paper by 
Castin [35J and works cited therein), but often with- 
out reference to the underlying supersymmetry of the 
Poschl- Teller potential. 

We can expand the deformation Sip over this or- 
thonormal basis set of eigenhmctions, 



5ip(x) = ao^>o(aO + ai^i(x) + / dka k ipk(x) 



(24) 



and then determine the coefhcients by projecting ( 18 1 
onto each eigenfunction: 



(E j + ^ aj = jdyr j (y)W(y). 



(25) 



Interestingly, the coefficient a± of the first excited bound 
state j = 1 remains undetermined because E\ = —1/2 
of Hq is exactly compensated by the +1/2. In other 
words, this mode appears as a zero-energy eigenmodc 
of the linear response kernel [T|[20 l f29 ] . This has a sim- 
ple, yet profound physical explanation: because ip\(x) — 
~yJi/Nd x ipo{x), the deformation ifo(x) + aiipi{x) = 
(fio(x + Sx) + 0(V 2 ) simply shifts the center of mass 
by 5x = —aiy/3/N. Thus we find, as argued in Sec. Ill 
above, that the soliton's center of mass is an indepen- 
dent dynamical variable that is influenced non-perturba- 
tively by the external potential. 

Let us then assume in the following that the soli- 
ton's center of mass has reached a position such that 
the external potential does not accelerate it anymore, 



namely that the right hand side of ( 25 ) for j = 1 van 



ishes. A sufficient condition for this is that the external 
potential V(z) is (locally) an even function around the 
soliton position z$. Then, the soliton shape deformation 
reads 



5ip(x) = a t/j {x) + / dka k tp k (x 



(26) 



with eq. ( 25 ) determining the coefficients ao and a k as 



linear functions of the external potential V(z), scaling 



by construction as V/|/x|. In Sec.[5]below, we will study 
the scaling with the wave vector for a simple lattice 
potential. 



4.3 Real-space response kernel 



By linearity, the soliton shape deformation ( 26 ) can also 



be conveniently expressed as a functional of the poten- 



tial (20) 



5tp(x) 



dyK(x,y)W(y), 



(27) 



with a symmetric kernel given by 



K(x, y) = -lUx)Mv) + 2 / dk M *WUv) , (28) 



The fc-integral can be evaluated, leading after some al- 
gebra to the following closed-form expression for the 
real-space response kernel: 



K(x,y) = - sech x tanh x sech y tanh y 



(29) 



[ cosh 2x + cosh 2y — | sinh 2x — sinh 2y| 
— 4 csch x csch y sinh \x — y\ — d>\x — y\] . 



Equations (27) and (29) together with (20) provide a 



pret-a-calculer expression of the soliton shape as func- 
tion of any given potential V(z). While soliton deforma- 
tion in external potentials is certainly not a new topic 
and has been considered in a rather large number of 
different contexts (see e.g. [5J and references therein), 
to our knowledge this useful expression of the linear- 
response kernel is new. 

This solution determines also the compressibility in- 
troduced in Sec. 4.1 From n(z) = [<j>o(z— zo)+S<j>(z)] 2 = 
n (z — z ) + 8n(z) + 0(V 2 ) it follows that 6n(z) = 
2<fio(z — zo)Sip(z), and we can read off from the previ- 
ous expressions that the density-density susceptibility 



of eq. ( 13 ) is given by 



x{z,z>) = 2ZMz-z*)K[ Z —^, Z —^)Uz'~z Q ). 



a ' s 



(30) 



The compressibility x{ z i z ') does not depend solely on 
z — z' , as it would in homogeneous systems. Instead, 
it refers to a distinguished point, namely the soliton 
position zo, and is only invariant under simultaneous 
translation of all 3 coordinates z, z', zq. 



4.4 Grand-canonical vs. canonical deformation 

The previous derivation used the grand-canonical set- 
ting at fixed chemical potential \x. The corresponding 
soliton deformation 5tp(x) =: Stp(x)\f 1 and density shift 
Sn(z) =: 6n(z)\u do not conserve the total number of 
particles. The change in particle number to order V/|/i| 
reads SN = J dzSn(z) = -// dzdz' X {z, z')V{z'). 

In order to calculate the soliton deformation at fixed 
number of particles N, i.e. to compensate SN, the chem- 
ical potential has to be adjusted. Since \fj,\ = g 2 N 2 /8, 
the required shift is S\fi\ — —2\fi\SN/N. Inserting p0| 
and performing the integral over z yields the relatively 
simple expression 



2| M | 



dy(l — y tanh y) sech(y) 



2 V(z 



<y) 



(31) 



Since the grand-canonical deformation Sip\n given by 



(27) is already of order V/ 1 ^ | , it is not affected by 



this shift to lowest order, but the original soliton back- 
ground is changed: (p (x)\ N = <A)(z)U + 9\fj.\(po(x)5\n\. 
Altogether, we find for the canonical soliton deforma- 
tion at fixed number of particles N 



2|// 



6ip(x)\ N = 5(p(x)\,i + 7-j — r(l - xtanhx)<po(x) (32) 



where Sip(x)\ fi is given by |27[ ), and (J|/x|/2|//| by ( |3T| ), 
with n = -g 2 N 2 /8, £ 



2/(N\g\) l mdMx) = V N / 2x 
secha; in all expressions. This formula, together with 

the kernel (29), constitutes the main achievement of 

this work. 



5 Special case: Lattice potential 

In order to illustrate the above results, we study in 
detail the case of a simple lattice potential with reduced 
wavevector fco = fc£: 



V(£x) = —Vo cos(k x) 



(33) 



with Vo > such that the center of mass of a soliton 
prepared at zq = (modulo the lattice period) sits in 
a potential minimum, where it remains classically. The 
limit fco <C 1 reduces to the particular case of a purely 
harmonic confinement studied recently by Castin [28 . 
Since the soliton deformation Sip is a linear functional 
of V, more general potentials (such as disordered ones) 
can be studied by applying the following results to their 
Fourier components. 




Fig. 2 (Color online) Relative shift in the chemical potential, 
5|/i|/|/i| as function of fco in a lattice of various depths Vo/ImI £ 
{0.2, 0.1, 0.05, 0.02, 0.01} [top to bottom]. Symbols: Numerical re- 
sults from an imaginary-time integration of the GP equation. 
Solid lines: linear-response prediction, eq. f''>4[ ). The inset shows 
the collapsed data, i.e. the raw data divided by Vo/ImI- A s fco>l 
or Vb/IH <S 1, the shift is correctly predicted and, as expected, 
small. 



5.1 Chemical potential shift 



First of all, the chemical potential shift (31) evaluates 

to 



%| Vp {TTk /2) 2 TTko 

\fx\ \n\ sinh(^ /2) 2 C ° S 2 ' 



(34) 



As fco = k£, — > 0, one finds S\/j,\ = Vo as expected, since 
this exactly compensates a global potential offset — Vq. 
Also not surprisingly, the shift vanishes as fco — > oo at 
fixed Vo since the condensate cannot follow these rapid 
variations. But it may perhaps come as a surprise that 



the behaviour as function of ko described by (34) is 
non-monotonic . 

We can compare these predictions to the results of a 
numerical integration of the imaginary-time GP equa- 
tion d T ip{x,r) = —H(p(x 1 r) that describes a steepest- 
descent trajectory towards the minimum of the canon- 
ical energy functional (i.e. (12) with /i = 0) [30,. This 
dynamics is not unitary, and the correct normalization 
to N must be re-established at each time step. When 
the stationary state is reached, tp(x, r) = e _/iT (^(x), the 
chemical potential is easily extracted from the required 
global renormalization factor. Figure [2] shows the mea- 
sured relative shift <5|//|/|/i as function of fco for various 
potential strengths (N = 100). The analytical predic- 
tion (34) is indeed found to be correct for Vo <C |/z|. The 
inset shows the collapsed data obtained after division by 
Vo/|/x|, clearly featuring the interesting fco-dependence, 
with the maximal shift reached close to fco = 1. 



5.2 Expansion coefficients and smaii parameter 



The expansion coefficients (25 1 read 



a 



&k 




f|(l + fc 2 )scch^, (35) 
~ 1 + fc 2 - 3fcg 

'Z\tT\ 4 (l + /c 2 )3/2(4+fc2)l/2 

x [sech | (A; + fe ) + sech §(fc - fc )] • (36) 

All coefficients scale as Vb/|/i| by construction. As func- 
tion of fcoj the bound-state coefficient ao ~ fcoe _7rfc °/ 2 
becomes exponentially small for a rapidly fluctuating 
potential &o 3> 1 (cf. the qualitatively similar behaviour 
of the effective potential ([5])) and thus does not con- 
tribute substantially in this limit. Away from |fc| = 
fco>l, also the coefficients a^ are exponentially small. 
However, the resonant amplitudes a±k ~ fc^" are only 
algebraically small. We therefore conclude that the modes 
±fco contribute dominantly, imprinting a pure sinusoidal 
wave onto the soliton, with an overall weight scaling as 



Vo 



2K, 
A; 2 






(37) 



In this smoothing regime, completely analogous to the 
repulsive case [21], the healing length (or chemical po- 
tential) drops out, and the small parameter of the ex- 
pansion rather is the ratio of potential amplitude to 
the kinetic energy Ef. associated with the spatial lat- 
tice wave vector k. 



5.3 Real-space response 



Unfortunately, even for the simple lattice potential ( 33 1 



the integral (27 1 over the kernel (29) cannot be evalu 



ated to simple closed form. However, it is easily calcu- 
lated numerically. Fig. [3] shows how the unperturbed 
soliton ground state <po(x) = \/N/2sechx with N — 
100 particles is deformed by a lattice potential with a 
rather large amplitude, Vo — 0.2\/i\ and intermediate 
lattice wave vector, kg — 5ir/8 ~ 1.9635. The predic- 
tion from the linear-response theory matches the full 
result quite well, considering that the perturbation is 
rather strong. To reach this agreement, it is essential to 
take into account the chemical potential shift accord- 
ing to eq. (32). The deviations d(f(x) are also shown 



in dashed. One clearly observes a concentration of den- 
sity around the potential minimum, as well as density 
depressions around the first potential maxima. 

To check the validity and precision of the linear re- 
sponse, we plot the deformation in direct comparison 
to the data from the numerical solution, together with 
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j \ — i fi(x)hn 

\ Sif(x) num _ 
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Fig. 3 (Color online) A soliton amplitude [green: ipo(x)] con- 
taining N = 100 particles and centered at x = is deformed 
by the lattice potential ( |33| > with Vo = 0.2|/j,| and wave vec- 
tor fco = 57r/8 [blue: V(a;)/|^|]. The numerically calculated GP 
ground state [black: ip(x) UUIn ] is very well approximated by the 
linear response result [red: </?(a:)iin] given by eq. ( 1326 . The respec- 
tive deviations 8<p(x) are also shown in dashed: the amplitude 
enhancement around the potential minimum at x = is clearly 
visible, just as the amplitude depression around the first potential 



maxima. 
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Fig. 4 (Color online) Deformation of the soliton shape by a 
lattice potential [blue: V(a;)/|/z|] of fixed wave vector fco = 
57r/8 and various amplitudes Vq. The linear-response result, eq. 
\'.V2\ [dashed red: <5ip(x)ii n ] approaches the numerically obtained 
ground state deformation [solid black: 8f{x) mlm \ as Vo/|H — > 0. 
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Fig. 5 (Color online) Deformation of the soliton shape by a lat- 
tice potential [blue: V(as)/|;t[] of fixed amplitude Vo = 0.1|/i| 
and various wave vectors fer> The linear-response result eq. 
(|32j [dashed red: Sip(x)n n ] approaches the numerically obtained 
ground state deformation [black: 8tp(x) nnln ] as fco 2> 1. 



the potential in Fig. H for fixed intermediate wave vec- 
tor &o = 57r/8 and various potential strengths, Vb/|/u| G 
{1,0.2,0.05,0.01}. By parity, we can restrict the plots 
to x > 0. As expected, the linear response approches 
the full solution very well as soon as V <C |/i|. 

Finally, Fig. [5] shows the same data for fixed po- 
tential strength Vq — 0.1|/i| and various wave vectors 
fco G 57r{jg, g, |, |}. Note that the vertical axis retains 
the same scaling in all four plots. As expected from 
the discussion in Sec. |5.2| above, the linear response be- 
comes increasingly accurate for larger k vectors, as the 
deformation itself becomes very small. The last set of 
data shows an almost pure sinusoidal deformation. In- 
deed, in the limit ko ^> 1, one can extract the leading 



contribution of (27) by partial integration and finds 



8ip(x) = — -ipo(x) coskox (38) 

Ek 

up to higher orders in Vo/Ek = 2Vo/k 2 . 



6 Concluding remarks 

In this article, we elaborate on two important aspects 
of bright cold-atom solitons in inhomogcncous external 
potentials. Firstly, we discuss the quantum dynamics 
of the soliton's center of mass in a disordered, spatially 
correlated potential. The center-of-mass wave function 
is predicted to show Anderson localization with a lo- 
calization length that can be estimated using the Born 



approximation of perturbation theory and is shown to 
be in an experimentally relevant range. 

Secondly, we calculate the imprint of a weak ex- 
ternal potential on the soliton mean-field shape in the 
ground state. For this, we determine the linear-response 
kernel, both for the wave function and the local den- 
sity, as well as the chemical potential shift. Amusingly, 
this can be done rather easily using the supersymmet- 
ric quantum mechanics of the Poschl- Teller or sech 2 - 
potential well. Finally, a detailed comparison to the nu- 
merically calculated ground-state solution is presented 
for a pure lattice potential with amplitude Vo and wave 
vector k. The small parameter of the linear response is 
Vo/|/x| in the regime fc£ ~ 1 of a potential varying on 
the scale of the soliton width. In the limit k£ — > 0, one 
recovers the case of harmonic confinement, studied re- 
cently by Castin [25] • Here, the small parameter for the 



soliton deformation (27) is Vbfc 2 £ 2 / 2 H = (w /2|m|) 2 in 



terms of the harmonic trapping frequency uj®. In the 
smoothing limit k£ 3> 1 of a rapidly fluctuating poten- 
tial, the condensate amplitude shows a pure sinusoidal 
imprint, the small parameter being Vo/Ek- 

Until now, we have discussed the center-of-mass quan- 
tum dynamics for fixed shape on the one hand, and the 
static properties of the mean-field ground state shape 
on the other. From the latter, we have learned that 
the impact on the soliton shape is indeed small for 
small Vo/|/l*|. In other words, the external potential can- 
not easily excite the internal modes since these have 
a gapped spectrum. Therefore, whenever the soliton 
moves slowly within a smooth potential, it is slightly 
polarized, but only adiabatically and reversibly (just 
as, say, an alkali atom is polarized in an optical dipolc 
potential). Therefore, we can expect the center-of-mass 
quantum dynamcis to be unharmed by shape excita- 
tions as long as the gap \/jl\ remains large compare to 
all other energies. 

It must be kept in mind, however, that Anderson 
localization is an interference effect relying on perfect 
phase coherence, at least over the time and distance 
needed to observe it. Already very small decoherence 
can kill this effect. For cold- atom solitons, the most 
dangerous sources of decoherence arguably are scatter- 
ing of background-gas atoms (cf. the analogous case of 
Fullerenes [31]) and three-body collisions, both of which 
can in principle be minimized in the experiment. How- 
ever, if the movement of a soliton inside a disordered 
potential were to radiate permanent excitations, as dis- 
cussed for the case of isolated point impurities in [5J, 
then this would constitute an intrinsic source of deco- 
herence. The detailed study of such effects, well beyond 
the simple mean-field estimate for the ground state pre- 
sented here, is left for future research. 
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